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ABSTRACT 

Spitzer IRS spectroscopy supports the interpretation that BP Piscium, a gas 
and dust enshrouded star residing at high Galactic latitude, is a first-ascent giant 
rather than a classical T Tauri star. Our analysis suggests that BP Piscium's 
spectral energy distribution can be modeled as a disk with a gap that is opened by 
a giant planet. Modeling the rich mid- infrared emission line spectrum indicates 
that the solid-state emitting grains orbiting BP Piscium are primarily composed 
of ~75 K crystalline, magnesium-rich olivine; ~75 K crystalline, magnesium-rich 
pyroxene; ~200 K amorphous, magnesium-rich pyroxene; and ~200 K annealed 
silica ('cristobalite'). These dust grains are all sub-micron sized. The giant planet 
and gap model also naturally explains the location and mineralogy of the small 
dust grains in the disk. Disk shocks that result from disk-planet interaction 
generate the highly crystalline dust which is subsequently blown out of the disk 
mid-plane and into the disk atmosphere. 

Subject headings: accretion disks — circumstellar matter — infrared: stars - 
planet-disk interactions — stars: individual (BP Piscium) 



1. Introduction 

BP Piscium (B P Psc) is an enigmatic h igh-Galactic latitude star enshrouded in gaseous 



and dusty material (jZuckerman et a/.ll2008l . hereafter Z08). Initially, BP Psc was thought 
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to be an isolated KO-type star undergoing its classical T Tauri star phase. Observations 
performed by Z08 instead found several indications that BP Psc is likely a first-ascent giant 
star that accretes material from a surrounding disk of gas and dust while launching bipolar 
jets. The first such sign that BP Psc might not be young was a lithium A6708 absorption line 
equivalent width that was ~7 times less than what would be expected for an early K-type 
classical T Tauri star. No known example exists in the literature of such strong lithium 
depletion for a Solar-like pre-main sequence star. A second indication that BP Psc was not 
young came from a low measured surface gravity. Although young stars typically have lower 
surface gravity than main sequence stars (as they are still contracting towards the main 
sequence), BP Psc's measured gravity of 10 2 ' 5 cm s -2 was significantly low er than those 
expec ted for 1-10 Myr old stars (which have gravity values >10 3 ' 5 cm s -2 ; see 



Baraffe et al. 



19981 ) and only served to further exacerbate the lithium situation. Yet another indication, 
although weaker than the first two, was a discrepancy between model and kinematically 
measured masses of BP Psc in the case that it was a classical T Tauri star. Such problems 
disappear if instead it is assumed that BP Psc is a first-ascent giant star. 

Although the sum of the observations performed by Z08 pointed towards BP Psc as an 
evolved rather than a young st ar, they were not nec essarily conclusive and further investiga- 
tion of BP Psc was warranted. iKastner et all ( 120081 ) performed molecular line observations of 
BP Psc showing that its molecular gas characteristics are compatible with t hose found in the 
expanding envelopes of yellow supergiant stars while IKastner et all (120 lOf ) imaged BP Psc 
with the Chandra X-ray Observatory fin ding a weak source with characteristics compatible 
with rapidly rotating G-type giant stars. iMelis et all (120091 ) reported the first unambiguous 
case of a dusty, accreting first-ascent giant star (TYC 4144 329 2), hinting that BP Psc 
could be part of a previously unidentified class of disk-bearing, accreting first-ascent giant 
stars (iMeha 120091 ) . Here we present Spitzer IRS spectroscopic measurements of BP Psc. 
These observations were designed to identify the constituents of BP Psc's dusty disk with 
the hope that such additional information would clarify the evolutionary status of BP Psc 
and, should BP Psc indeed be a first-ascent giant star, provide clues as to how it may have 
been rebirthed into characteristics typically associated with young, planet-forming stars. 



Observations 



BP Psc was o b served with the S pitzer Space Telescope Infrared-Spectrograph (IRS; 
Werner et all 12004 ; iHouck et all 120041 ) in stare mode on 2007 August 02 with the low- 
resolution modules affording resolving powers of ~60-120. Since BP Psc is such an infrared 
bright target, 1 cycle of 6 seconds ramp duration for each of the two nods was sufficient 
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to obtain well-exposed (S/N>100 per pixel) spectra for each of the Short-Low (SL) and 
Long-Low (LL) modules. Data were pre-processed with Spitzer Science Center IRS Pipeline 
Version S16 . 1 . for the SL modules and S17 . 2 . for the LL modules. Subsequent bad pixel 
masking and interpolation was performed using IRSCLEAN VI. 9 and a "master" bad pixel 
map. The master bad pixel map was generated by running the interactive IRSCLEAN_MASK 
code on each of the four unique basic calibrated data spectral images per module, then 
combining these four output files and the campaign rogue mask that identifies known bad 
pixels. The purpose of the master bad pixel map method was to identify and remove bad 
pixels that might otherwise go unnoticed under the well-exposed spectrum. 

Individual nods from a pair were subtracted, then these background-subtracted spec- 
tra were extracted and flux-calibrated using the SPICE software package. Each extracted, 
calibrated nod spectrum was compared against its nod pair to ensure fidelity in the final 
spectrum. Data from SL order 3 (the so-called "bonus order" ) were removed due to a poor 
match between the two nods. Nod-averaged spectral orders were generally well-matched in 
terms of flux-continuity. Minor adjustments, where necessary, were made such that an indi- 
vidual order matched the flux of the bluer order. The final spectrum is displayed in detail in 
Figure [T] and is overplotted on a complete spectral energy distribution for BP Psc in Figure 

m 

On 2007 July 08 BP Psc wa s serendipitously imaged at 24 /xm with MIPS (Multiband 



Imaging Photometer for Spitzer; iRieke et al\ 12004 ) in a program to monitor the dust cloud 



following Earth's orbit (AOR 18701312, PI Jayaraman; BP Psc did not fall into the 70 fim 
field of this program). We performed aperture photometry on the 24 ^m post-basic calibrated 
data mosaic image produced by the Spitzer Science Center MIPS pipeline (version S16 .1.0). 
We used the Spitzer Science Center recommended aperture correction of 1.167 for an aperture 
radius of 13" and sky annulus inner and outer radii of 20" and 32 ", respectively. We detect 
a flux density of 5.7 Jy at 24 /iin with an uncertainty of 10%, consistent with the IRAS 25 
fim detected flux and the IRS flux (see Figure |2]). 



3. IRS Modeling 

Before discussing in detail IRS spectroscopy modeling and results it is necessary to recall 
what is known about BP Psc's disk (see Z08 for details). Resolved molecular gas emission 
detected through Submillimeter Array aperture synthesis imaging indicates that BP Psc 
has a gas disk in Keplerian orbit. Keck II adaptive optics imaging resolved a dust disk in 
scattered light. Near- to far-infrared photometric measurements of BP Psc reveal that ~75% 
of BP Psc's observed light is reprocessed into the infrared and that the orbiting material 
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is well modeled as residing in two regions having temperatures of ~1500 and ~210 K (see 
Figure [2]). From a Submillimeter Array continuum detection, Z08 calculate that 7 M e of 
dust orbits BP Psc in the case that it is a giant star 300 pc distant from Earth. Synthesis 
of the imaging data led Z08 to suggest that we view BP Psc's disk close to edge-on (i~75°) 
and that the disk has a flared morphology. 



3.1. Disk Structure 



Acke et al\ (120091 ) and iFurlan et all (120061 ) validated the use of IRS-based spectral in- 



dices to study disk structure by comparing indices computed from data to those computed 
from disk models. We seek to gain insight into BP Psc's disk str ucture by comparin g its 
IRS spectral indice s to those computed for the disks studied by Acke et al\ (120091 ) and 
Furlan et al\ (120061 ) . We calculate the value of the following indices: an excess at 7 jum 
of 4.55 magnitudes, the ratio of flux at 13.5 /mi to that at 7 /mi ([13.5/7]) of 0.69 magni- 
tudes, and the ra tio of flux at 30 /im to that at 13.5 /im ([30/13.5]) of 0.96 magnitudes (see 
Acke et a/.ll2009i . for how these values are derived); the ratio of integrated continuum flux 
near 6 /im to that near 13 /tm (p.Q-13) of 2.10, the ratio of integrated continuum flux near 
13 /tm to that near 25 /tm (n 13 _ 25 ) of 1.93, an d the ratio of integr ated continuum flux near 



6 /tm to that near 25 /im (n 6 _ 2 s) of 2.03 (see IFurlan et a/.ll2006l . for how these values are 
derived). These values suggest that BP Psc has strong near- infrared excess emission and 
steeply rising flux towards longer wavelengths (consistent with BP Psc's SED; see Figure [2] 
herein and Figure 6 of Z08). As mentioned above, the real power of these indices to con- 
strain the structure of BP Psc's disk comes from comparing them to indices derived for other 
well-modeled disk-bearing objects. Such comparisons are presented in Section [U 



3.2. Dust Composition 



The IRS spectrum of BP Psc shows several broad emission lines from solid state tran- 
sitions (Figures [U|3]). To model these features, and hence deduce the nature of the emitting 
particles, we proceed in fitting the spectrum via two paths. In one path we attempt to iden- 
tify grain species by subtracting the most obvious contributor and searching for residuals 
indic ative of additional grain species (as was done successfully for HD 100546 by lMalfait et al. 
19981 ). During the se subtraction fit s we try various grain shapes and sizes (using the absorp- 
tion coefficients of iMin et all 120071 ) to identify the best shape distribution and grain size to 
use. In para llel with these subtr action fits, we fit the data using the x 2 minimization routine 
described in lGielen et al\ (120081 ). We refer the interested reader to their work, and references 
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cited therein, for details about model inputs. The strategy employed by iGielen et al\ (120081 ) 
assumes that solid state emission features originate in an optically thin disk atmosphere. 
They attempt to fit mid-infrared spectra with a blackbody model for the disk continuum 
emission and an emission feature model computed from the linear combination of dust ab- 
sorption coefficients multiplied by different blackbody temperatures. The x 2 of the model 
fit to the data is then minimized by adjusting various parameters; e.g., grain temperature, 
grain species, and intensity of grain emission. Both the preliminary x 2 fit an d subtraction 
fits indicate that the grains producing the emission features in BP Psc's spectrum are sub- 
micron (~0. 1-1.0 jttm) in size. We also find that both m ethods obtain consi stent results in 



grain species identified. Similar to what was found by IGielen et al\ ( 120081 ). we find that 



dust emission features cannot be reproduced with spherical grains (e.g., Mie theory), but are 



(Shkuratov & Grvnkol 


2005 


Min et al. 


e Approximation ( 


Drainel 


19881). It is 



noted that the grain shape approximations used sometimes do not reproduce exactly the 
grain emission featur es (and hence likely are not a complete match to the true grain shapes); 
see the discussion in IGielen et al\ ( 120081 ). In the following three paragraphs we describe the 
resulting grain species identified from these two methods. 

Three features — with wavelengths near 33, 28, and 23 /xm — stand out in BP Psc's IRS 
spectrum. These features are typically associated with crystalline, magnesium-rich olivine 
grains (forsterite, Mg 2 Si04). Fits to the forsterite features left significant residual emission 
in the 20-35 //m range. Experimentation with various grain species indicates a contribution 
from crystalline, magnesium-rich pyroxene (hereafter enstatite, MgSiOa). Fits including 
enstatite left little residual in the 23-35 /im range. 

Our forsterite and enstatite fits alone could not reproduce the feature near 10 /im. 
The ~10 /im feature is typically associated with amorphous dust of the pyroxene or olivine 
species. Experimentation with the two suggested that amorphous, magnesium-rich pyroxene 
(hereafter pyroxene, MgSiOa) provided a better fit. 

The remaining residual emission — near wavelengths of 9, 13, 16, and 21 /im — is 
reminisce nt of the features pr oduced by annealed, crystalline silica (hereafter cristobalite, 
SiO 



2', see 



Sargent et a/.ll2009bl . and references therein). It is noted that the 9 and 16 micron 
complexes are marginally detected in the residual spectrum. In light of this fact, and in 
the interest of completeness, we reviewed other polymorphs of silica to see if they could 
reproduce the remaining spectral features, a- quartz is ruled out by it s characteristic 25 /im 
feature which would be quite strong if present (IGervais fc Pirioull 19751 ). /3- quartz is ruled out 



by the high temperatures necessary to excite its resonances (Td us t>850 K; 



Gervais & Piriou 
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19751 ); our above analysis suggests the small dust grains seen in BP Psc's IRS spectrum are 
significantly cooler (see also below). The polymorphs coesite an d stishovite are also ruled out 
due to their unique emission spectra (see lSargent et a/.ll2009bl and references therein). This 
leaves cristobalite, amorphous silica, and obsidian as potential silica polymorphs to explain 
the residual emission. Since obsid ian and amorphous silica have similar spectral features (see 
Figure 3 of ISargent et all l2009bl ). we use amorphous silica to repres ent both p o lymor phs. 
To model silica emission features we adopt optical constants from iMin et all ( 120071 ) for 
amorphous silica (hereafter silica) and from C. Lisse ( 2009, private com munication; see the 
discussion in the Supplementary On-Line Material to iLisse et all 120061 . for a description of 
the input optical constants) for cristobalite. Absorption coefficients are calculated for silica 
using GRF and for cristobalite using a continuous distribution of ellipsoids (CDE). Models 
using silica result in higher x 2 than do fits using cristobalite and fits without any SiC>2 
polymorph provide yet even higher x 2 values; a comparison of fits using cristobalite, silica, 
and neither of cristobalite or silica are shown in Figure HI As such, we consider SiC>2 to be 
significantly detected in the IRS spectrum of BP Psc and conclude that annealed silica is 
the most likely polymorph. 

Some of the above mentioned transitions are coincident with those of polycyclic aromatic 
hydrocarbons (P AHs). PAHs, however, typically display strong 8.0 //m features (see e.g., 
Sloan et a/.ll2007l . and references therein). No such feature is seen in our spectrum of BP Psc, 
thus we conclude that PAHs do not significantly contribute to the strength of the observed 
solid state features. 

We briefly mention the affect of reddening on mid-infrared solid-state emissio n features. 
Highly reddened stars can exhibit depressed 10 /mi features (IGielen et a/.ll2008l ). BP Psc, 
at a galacitc latitude of —57.2°, is unlikely to be experiencing any significant interstellar 
extinction regardless of its distance. The dusty material in orbit around this star could in 
theory provide substantial reddening if the grains are small. In their analysis, Z08 compared 
low-resolution optical spectra of BP Psc to non-dusty stars of similar spectral types and 
showed that BP Psc is not reddened along our line of sight. Thus, reddening is unlikely to 
contribute to the observed strength of the 10 /mi structure. 

With rough grain parameters identified we proceed in a final x 2 minimization fit to the 
spectrum. We restrict input dust species, grain shapes, and grain sizes to those as determined 
from the subtraction and initial x 2 fits. We allow as free parameters three independent 
dust continuum temperatures, three independent dust species temperatures, two dust grain 
sizes (~0. 1-1.0 /an and ~2. 0-4.0 //m), and five dust species — forsterite, enstatite, pyroxene, 
amorphous olivine, and cristobalite. W e note that the final x 2 minimization routine employed 
herein differs slightly from that used by lGielen et all (120081 ). Our routine assigns grain species 
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into one of three distinct groups where each group is allowed its own dust temperature and 
temperature step size in the model grid. Amorphous materials (pyroxene and olivine) form 
one group which has a temperature grid step of 25 K, crystalline olivine and pyroxene form 
another group which has a grid step of 25 K, and silica is in its own group that has step 
size of 50 K. The final model fit is displayed in Figure [3j Temperatures, mass fractions, 
associated errors, and an estimated mass per spe c ies ar e reported in Table [TJ Errors on 



these parameters are calculated as in lGielen et all ( 120081 ): we randomly add gaussian noise 



with a sigma distribution determined from the spectrum statistical uncertainty per pixel at 
each wavelength point. In this manner we generate 100 synthetic spectra, all consistent with 
our data, on which we perform the same fitting procedure. The resulting slightly different 
fit parameters are used to derive the mean (our best-fit value) and the standard deviations. 
The uncertainty on the mass absorption coefficients are not taken into account in this error 
determination. 

Here we summarize some of the main results of the spectral fitting. No larger sized 
(~2. 0-4.0 ^m) grains are required in the fit. Amorphous olivine (Mg2Si04, hereafter olivine) 
is not necessary to obtain the minimum x 2 model fit. SiC>2, however, is required despite its 
low mass fraction and its inclusion in the fit results in a sign ificantly reduced final x 2 (~84 if 



included and ~110 if not; see Section 7 of iGielen et a/.ll2008l . for how we calculate x •) From 



Tabled] it is clear that the crystallinity fraction of the small dust grains is high, ~96%, with 
the majority of this coming from forsterite and enstatite. 

The sub-micron sized dust grains orbiting BP Psc, in the case that it is a first-ascent 
giant star, would be radiatively blown-out without a spatially coincident gas disk (this would 
not be the case if BP Psc were a T Tauri star). To illustrate this we assume BP Ps c has 



a luminosity of 100 L Q and a mass of 1.8 M Q . With Eq. 1 from IChen et all (120061 ). and 
assuming grain densities of 3.3 and 2.3 g cm -3 for silicates and SiC>2 respectively, we estimate 
that grains smaller than ~10 /xm in size will be radiatively blown-out of the BP Psc system 
in the absence of gaseous material. 



3.2.1. Unmodeled Features 



Identification of the ~6.5 fim absorption/emission complex in BP Psc's IRS spectrum 
(see Figure [T]) is complicated by its weak signatu re. It is possible that this feature arises 



from the v<i vibrational band of water (see e.g., IWoodward et al\ 120071 ). This particular 
fundamental of gas-phase water presents rovibrational emission features near 6.5 /im similar 
in morphology to the weak features obser ved at that location in BP Psc's IRS spectrum 



(compare Figure [T] herein and Figure 2 of IWoodward et all 120071 ). Since this is not a high 
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fidelity identification, we will refrain from attempting to interpret the presence of water 
vapor in BP Psc's disk. High-spectral resolution, high signal-to-noise spectra around 6 /im 
will be necessary to reliably identify and characterize this feature. 

There is potentially unmodeled solid-state emission remaining near 20 /im. If the 20 
jum emitter is a dust species not already present in the model fit, then it would be unlikely 
to have other transitions in the ~5-40 /im range as the rest of BP Psc's IRS spectrum 
is reproduced well by the model di splayed in Figure pH A potential carrier of the 20 /im 
emission feature is FeO nano dust ( jZhang et a/.ll2009l . and references therein). FeO nano 
dust hosts one solid-state transition near 20 /im. Temperature controls the width and peak 
wavelength of the FeO nano dust feature, with temperatur es of 5,100 K shifti ng the feature's 



peak wavelength to 20.1 /im and its FWHM to 2.4 /im (IPosch et all 12004. and ref erences 
therein). Such a peak position and profile (see also Figure 9 of IPosch et all 12004 ) could 



potentially account for the unmodeled emission in BP Psc's IRS spectrum. FeO nano dust 
can only exist in environments with low ultraviolet photon density and cool temperatures; 
ultraviolet photon s will dissociate F eO while high temperatures will prevent Fe and O from 
binding together ( IPosch et all 12004 ). Fe and O orbiting BP Psc at similar locations as the 
other dust species would be sufficiently cool to bind together while the late spectral type of 
BP Psc could provide the requisite low ultraviolet photon density to allow the FeO to stay 
bound. 



4. Comparison to Other Circumstellar Disks 

In the following subsections we compare BP Psc's mid-infrared spectrum to those of 
other well studied disk-bearing objects in an effort to better understand the physical envi- 
ronment responsible for BP Psc's disk structure and composition. 



4.1. T Tauri Stars 



T Tauri stars are pre-main sequence, Solar-like stars that are still surrounded by sub- 
stantial amounts of primordial gaseous and dusty material. The IRS-indices (see Section 
3.ip for these disk s are consistent with flat, optically thick, active accretion disks (see e.g., 
Furlan et a/.ll2006l). BP Psc's IRS-ind ices are strongly di s crepa nt with those for the T Tauri 
stars presented in lFurlan et all (120061 ) and lWatson et all (120091 ) — where the Taurus median 
of n 6 _ 13 , n 13 _25, and n 6 _ 2 5 are —0.82, —0.17, and —0.48 respectively — suggesting that the 
structure of BP Psc's disk is significantly different from T Tauri star disk structures (see 
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Section 

T Tauri star dusty material shows a variety of compositional characteristics that are at 
odds or incompatible with what is observed for BP Psc: 



• T Tauri stars sometimes host disks with high crystallinity fractions, but these grains are 
typically accompanied by significant quantities of a morphous species and are confined 
to within <10 AU separation from their host stars (IWatson et a/.ll2009l ). Reasonable 
age and mass estimates of BP Psc as a classical T Tauri star (see Table 6 of Z08) would 
require the cool forsterite and enstatite grains to reside at separations >12 AU. 

• If BP Psc were a young T Tauri star, its disk atmosphere would have one of the highest 
masses in small grains (~0.1 Lunar masses if located 80 pc from Earth; see notes to 
Table 1) for a star of its class (e.g., ISargent et al\ l2009al Jbl. quote a range of mass in 
small dust grains for the T Tauri stars they studied of ~0.8-700xl0 -4 Lunar masses). 



Watson et all (120091 ) find in the T Tauri stars they studied an anti-correlation between 
the strength of the 11 /mi crystalline olivine featur e and the tota l stren gth of the entire 
10 jum solid-state emission complex. Additionally, IWatson et all ( 120091 ) identify a posi- 
tive correlation between the total strength of the 10 /im complex and the magnitude of 
the ni3_25 index. These correlations suggest that stronger crystalline emission features, 
and hence more chemically evo l ved gra in populations, accompany disk settling in T 
Tauri star disks. ISargent et all (j2009al ) found that chemical evolution doesn't neces- 
sarily imply grain growth h as occurred, further suggestin g that chemical evolution and 



sedimentation are tied (but iDullemond fc Dominik 



2008. show that dust sedimentation 



alone, without grain growth, cannot reproduce the observed correlations). BP Psc's 
grains, despite their nearly 100% crystallinity fraction, are still suspended in the flared 
disk atmosphere, apparently in contrast with the observed T Tauri star trends. 



4.2. Herbig Ae/Be Stars 

Herbig Ae/Be (HAeBe) stars are the intermediate mass counterparts of classical T Tauri 
stars. HAeBe stars are much hotter and more luminous (median T n ff~6000 -7000 K, L~10- 
1000 L Q depending on their progress along pre-main sequence tracks; see iPalla fc Stahler 



19931 ) than lower mass, solar-like T Tauri stars and hence provide a far different circumstellar 



environment. 



BP Psc 's spectral energy d istribution is reminiscent of Meeus group I HAeBe stars (see 
Figure 4 in iMeeus et a/.ll200ll ). The IRS-indices for BP Psc, when plotted in Figure 2 of 
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Acke et all (120091 ) , lie in the region of infrared excess space describing HAeBe disks having a 
Meeus group II disk structure (where group II sources have an outer disk which is protected 
again st direct stellar radiation by the puffed-up inner disk rim; see iDullemond fc Dominik 
2004 and references therein). This behavior is peculiar, but repeated in a well-studied 
HAeBe star, HP 1 00546. This source is given as an exemplar Meeus group I source in 
Meeus et all (120011 ). However, it too resides in the region of Figure 2 of lAcke et all (12 009) 



descri bing Meeus group II disk structures. Recent interferometric work by iBenisty et al. 



(120101 ) has confirmed that HD 100546 has a gap in its disk (see Section 16. 2p . It would 
seem that gapped disks can provide confusing results when relying on spectral indices alone. 
Nonetheless, BP Psc's spectral energy distribution, and hence disk structure, is consistent 
with those found for HAeBe stars. 

The match of BP Psc's disk structure with those of HAeBe stars suggests that the 
inner disk morphology for both types of disk systems is determined from the same physical 
effect. It is not yet clear what determines this inner radius environment, but the ubiquity of 
hot excess temperatures of ^1500 K that fit the near- infrared "bumps" in HAeBe spectral 
energy distributions (thought to arise from puffed-up inner disk rims) suggests that dust 
sublimation is responsible. This is in contrast to Solar mass T Tauri stars, whose inner disk 
morphologies are determined by magnetosphe ric accretion (which act on T Tauri disks at 
~7 R*, close to sublimation temperatures; see lBouvier et all 120071 . and references therein). 



ISO, IRS, and ground-based mid-infrared spectroscopy of HAeBe stars indicates disks 

sometimes populated by PAHs and/or significantly evolved dust as evidenced by high mass 

fractions of larger than sub-micron s i ze, chemically and thermally processed gra ins (jBouwman et al. 
200 ll ; lAcke k van den Anckerlbooi Ivan Boekel et a/.ll2005l : iKeller et a/1l2008[ ). One HAeBe 
star, HD 100546, shows dust emission features in its mid- infrared spectrum that are sim- 
ilar to those seen i n our IRS spectrum of BP Psc (Figure [2]). HD 100546 was shown by 
Malfait et all (119981 ) to have strong contributions to its mid-infrared spectrum from forsterite 
and PAH emission. The grains or biting HD 100546 ar e sub-micron sized and originate from 
two distinct temperature regions (IMalfait et a/.lll998l ). More discussion of HD 100546, and 
its relation to BP Psc, can be found in Section [61 



4.3. First-Ascent Giant Stars 



Three other first-ascent giant stars have well characterized i nfrared excess emission. Two 



of these objects, HD 233517 (a K2 III star; iFekel et al 
and HD 100764 (a ~4400 K, R-type Carbon star; Iskinner 



996 



Jura 



1994 



2003 



Sloan et al 



published IRS spectra. The third, TYC 4144 329 2 (an F2-type star; Melis et all 120091 ). has 



ura et a/.ll2006f ) 



2007]) , also have 
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existing IRS spectroscopy that will be mentioned here (a more detailed description of TYC 
4144 329 2's IRS spectrum will appear in C. Melis et al. 2010 in preparation). The infrared 
excess emission for these three objects is modeled as arising from a Keplerian disk of dusty 
material that orbits the central giant star. In the case of HD 100764 and T YC 4144 329 2 , 
near-infrared exc ess emission is present indicating dust as hot as ~1500 K ( jSkinnerl 1 19941 ; 



Melis et a/.ll2009l ). HD 233517 exhibits infrared excess emission indicati ng a flared , cool dust 



disk (T dust ~70 K) orbiting at ~45 AU separation from the giant star ( 1 Jural 120031 ). 

IRS spectroscopy of all three of the above mentioned giant stars indicates that they are 



orbit ed by PAHs and small, amorphous, silicate dust grains (IJura et a/.ll2006l ; ISloan et al 



20071 . C. Melis et al. 2010 in preparation). It is interesting to note that all three of these 



giant stars host obvious PAH features while BP Psc does not (see Figure d} although the 
~6.5 fxm feature in BP Psc's IRS spectrum remains unidentified, it is unlikely to be the 
characteristic 6.2 /iin PAH transition). 



4.4. Post-AGB Binaries 



With the IRS and TIMMI2 mid-infrared spectrographs, iGielen et al\ (120081 ) performed 
a comprehensive study of the circumstellar material in Keplerian orbits around post-AGB 
binaries having a median stellar effective temperature of ^6000 K and assumed luminosities 
in the range of ~5000±2000 L Q . T he disks around p ost-AGB binary systems are well- 
modeled as passive, irra diated disks ( Deroo et al. 2007al lbl) . Some of the post-AGB binary 
IRS spectra presented in lGielen et al. ( 2008 ) have similar slopes as the BP Psc IRS spectrum, 
indicating that they could have similar structure. 



Gielen et all (120081 ) found that the stable disks orbiting post-AGB binaries host popula- 



tions of hot and cold amorphous and crystalline grains. Grains in post-AGB binary dust disks 
show evidence for strong processing, both in the form of grain growth beyond sub-micron 
sizes and crystallization. The physical process by which disks form in post-AGB binary 
systems is not well understood, but it is e xpected that intera ctions between the evolved star 
and its companion play a necessary role (IDeroo et q/.ll2007bl ). 



Similar to our discussion of radiative blow-out for grains orbiting BP Psc, it can be shown 
that post-AGB binaries must also have significant quantities of gas spatially coincident with 
their dusty material. Since the minimum grain size that can orbit a star scales linearly with 
luminosity, and since post-AGB binary systems have similar masses as first-ascent giants, 
we estimate that the minimum size grains that can orbit a post-AGB binary in the absence 
of gas would be >100 fxm in size. 



-12- 



5. Implications for the Evolutionary State of the BP Psc System 



Z08 established through various methods that BP Psc has the effective temperature of 
a G8-K0 star and likely has the luminosity of a first-ascent red giant star. We seek to test 
the luminosity result under the assumption that the effective temperature is correct. 

Comparison of BP Psc's dust mineralogy to those of other classes of disk-bearing objects 
does not provide any significant insight into BP Psc's evolutionary status. Although there 
are discrepancies with what is observed for classical T Tauri stars, there are also discrepancies 
with the constituents observed in disks around first-ascent giant stars; for example, other 
giant stars exhibit detectable PAH emission whereas BP Psc does not. BP Psc, with an 
accretion rate of 10~ 8 M yr -1 (Z08), almost certainly emits sufficient ultraviolet radiation 
to excite PAH emission; especially if late-type giant stars like HD 233517 and HD 100764, 
who do not have any obvious accretion, can excite PAH grains. The absence of obvious 
PAH emission suggests that BP Psc lacks abundant PAH molecules (perhaps the effect of 
carbon-poor disk material) and/or that the strong near-in frared continuum emission from 
BP Psc's inner disk rim veils PAH emission features (e.g., Geers et al. 20061 ). In terms of 
the presence of a crystalline dust population and lack of PAHs, BP Psc's dust characteristics 
most resemble those of post-AGB binaries (it is noted that Z08 demonstrated that BP Psc 
cannot have a lumi nosity as large as t hat of an AGB star). BP Psc's dust differs from the 
objects studied by iGielen et al\ (120081 ) as the majority of BP Psc's grains are sub-micron 
sized and cool (Td us t<300 K) while the majority of post-AGB binary stars exhibit hot grains 
(Tdust^300 K which results in strong emission features in the ~10 /im region) with sizes 
larger than ~1 /im. 



As discussed in Sections [3] and \A.2\ BP Psc's disk has striking structural similarities to 
disks around HAeBe stars; in particular, the necessity for the puffed-up inner disk rim that 
is characteristic of HAeBe disks. This result suggests that BP Psc, if a classical T Tauri 
star, lacks strong magnetic fields that effectively funnel material onto the star and prevent 
the build-up of a substantial inner disk rim. It is hard to imagine BP Psc as a classical T 
Tauri star without strong magnetic activity (e.g., iPreibisch fc Feigelsonl 120051 ). Measuring 
Zeeman splitting of BP Psc's atmospheric absorption lines could provide concrete evidence 
against t he system being; pr e-main sequence should the lack of a strong magnetic field be 
verified (IKastner et a/.ll2010l . favor a scenario where BP Psc's X-ray emission is generated 
through either coronal processes or star-disk interactions, both of which require BP Psc to 
be magnetically active). 
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6. BP Psc's Disk Origin and Evolution 



Analysis of our IRS data suggests that BP Psc is more likely an evolved sta r than a T 
Tauri star, consistent with the results obtained by Z08 and iKastner et all (120081 ). How did 
this giant star come to be orbited by so much dusty and gaseous material? Z08 postulate 
that BP Psc, if an evolved star, could have been rebirthed into its young star-like dusty 
and gaseous environment by consuming a short orbital-period companion at the onset of 
its post main-sequence evolution. Such a d isk formation me chanism is also fav ored for th e 
dusty first-ascent giants TYC 4144 329 2 jMelis et a^.l l2009h and HD 233517 J Jural l2003h . 
Unfortunately, the literature is lacking binary engulfment models that produce results con- 
sistent with observed outcomes — giant stars bearing disks in excess of their stellar radius 
and no evidence for the hypothetical consumed companion. The model that provi des the 
most reasonable match to the final observed outcome is outlined in lShu et al\ (I1979I ). Their 
analytical model entails interaction between an evolving star and its short-orbital period, 
less massive secondary that causes matter to leave the system through the L2 lagrangian 
point in a spiral pattern. Eventually the spiral outflow of material will return to strike itself 
resulting in the formation of a stable disk of material that accretes onto the central giant 
star and whose maximum size is governed by conservation of angular momentum. 



6.1. Does BP Psc's Dust Composition Result from Binary Engulfment? 



According to the IShu et all (119791 ) model, material lost from binary interaction will 
initially be close to the host giant star and thus would be subjected to the intense giant 
star radiation field. Once this material passes beyond the sublimation radius it can begin 
to coalesce into dust grains. These dust grains should be efficiently annealed with the heat 
from the nearby luminous star, resulting in a high fr action of dust be i ng cry stalline. In their 
study of the 10 jum emission feature of HAeBe stars, iBouwman et al\ (120011 ) found that with 
increasing abundances of forsterite one expects to see larger mass fractions of crystalline Si02- 
Such a result is indicative of the thermal annealing process which transforms amorphous into 
crystalline grains with Si02 as a by-product. That we see Si02 in the material s urrounding 



BP Ps c lends some support to a disk formation model similar to that presented by lShu et al. 



(119791 ). The orbital semi- major axis of this crystalline material will continue to expand until 
the disk reaches its maximum size. Thus, highly crystalline material would be transported 
to the cooler, outer regions of the disk where it is observed. 

There are some problems with the above grain evolution scenario. One is that other first- 
ascent giants with dusty material thought to originate in a manner similar to BP Psc's do 
not show evidence for highly crystalline material. A second problem is that the polymorph of 
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Si02 most likely present in BP Psc's IRS spectrum, cristobalite, must be cooled very quickly 
after being annealed or it will revert to lo wer temperature polym orphs that do not exhibit 
the characteristic 16 //m emission feature (jSargent et a/.ll2009bl ). It is difficult to imagine 
that the disk excretion could be slow enough to efficiently anneal Si02 into a cristobalite 
polymorph and at the same time rapid enough to transport this cristobalite to cooler regions 
before it can revert. More detailed modeling and simulations of disk formation through 
engulfment could further explore these issues. 



6.2. A Planetary Model to Account for BP Psc's Observed Disk Structure 

and Composition 

It seems unlikely that the composition of the dust orbiting BP Psc is linked to its disk 
formation mechanism due to the problems discussed above. Instead, the dust mineralogy 
could be the result of evolutionary processes which occur in protoplanetary disks and are 
understood to be the precursors to (or indicators of) planet /planetesimal formation. Under 
such an assumption our IRS results could provide insight into whether or not new planetesi- 
mal formation is occurring in BP Psc's disk. As discussed above, the presence of cristobalite 
would require an event that heated the disk to temperatures of ~1000 K and then allowed it 
to cool rapidly. Such an event likely took place in BP Psc's outer disk where th e disk equilib 



rium t emperature is suf ficiently low to preven t the cristobalite from reverting (jSargent et al. 



2009bl ). As discussed in lSargent et all ( ]2009bl ) . disk- lightning or nebular shock scenarios can 
explain the presence of cool, highly crystalline material (especially cristobalite). However, 
they do not explain why so much of this material would be in BP Psc's flare d disk atmo- 



sphere (lightning and shock s are expected to mainly act in the disk interior; iPilipp et al. 



19981 ; lHarker fc Deschll2002l ). especially when one considers that grain processing shoul d be 



linked with grain growth and sedimentation (e.g. 
references therein). 



Natta et al. 2007: Watson et al. 2009, and 



Solutions to these problems come from comparing BP Psc to HD 100546. iBouwman et al. 



( 120031 ) suggest that HD 100546, who presents a spectral energy distribution (Figure [2]) and 
crystalline silicate emission similar to those of BP Psc (with the exception of the presence 
of SiC^), has a proto-gas giant planet carving a gap within its disk. Such a gap naturally 
explains the enhanced mid- to far-infrared excess emission relative to the near-infrared excess 
emission in HD 100546 when compared to other HAeBe stars as the far side of the gap would 
be directly irradiated by the star. This direct stellar irradiation would produce a vertically 
extended wall in the disk, allowing it to intercept a large fraction of stellar light which would 
then be reradiated as blackbody emission having the characteristic temperature for the disk 



- 15 - 



gap location (~200 K for HP 10054 6 and BP Psc; see Figure [2]). Interferometric observa- 
tions of HD 100546 by lBenisty et all (120101 ) confirm that HD 100546 does indeed have a gap 
within its disk. Recent optical and ultraviolet spectroscopic studies of HD 100546 further 
supp ort the interpretat i on that this star has a giant planet orbiting within the gap in its 
disk jGradv et al.ll2005[ kcke fc van den Anckerl l2006h . Since BP Psc and HD 100546 have 
similar SEDs (Figure |2]) and hence similar disk structure, it is reasonable to surmise that 
the same mechanism responsible for HD 100546's disk structure operates also in BP Psc's 
disk. Thus, it is plausible that a giant planet orbits and opens up a gap within BP Psc's 
disk. Such a giant planet would have an orbital semi-major axis of ~4 AU (assuming BP Psc 
has a luminosity of ~100 L Q ). We note that theoretical works predict that any such object 
have a mass that is <10 Mj up as more massive bodies would restri ct the flow of material 
towards the star resulting in a depleted inner disk (ILubow et a/.lll999l ). Such depletion is not 
compatible with observations that indicate large quantities of gaseous and dusty material in 
BP Psc's inner disk (Z08). It is not clear whether this giant planet would be from BP Psc's 
original planetary system or if it would be a protoplanet forming from the material in BP 
Psc's giant star disk. 

Two mechanisms to explain the observed dust mineralogy are discussed in the following 
subsections. Both rely on the existence of a giant planet opening a gap within BP Psc's disk. 



6.2.1. Colliding Planetesimals 



The giant planet orbiting within BP Psc's disk could stir up rocky bodies that also 
orbit in the disk gap. These rocky bodies could be remnant planetary building material 
left over from BP Psc's first planetary system that orbit at similar stellar separations as 
the giant planet (like Jupiter's trojan asteroids) or they cou ld have formed in-situ in BP 
Psc's revived giant star disk. Following iBouwman et al\ (120031 ). we expect ensuing collisions 
between these stirred planetesimals to release previously generated crystalline material. This 
collisional residue is then blown out of the disk gap by various processes (see Section 16.2.21) 
and eventually lands in the disk atmosphere where it is observed. 

If the colliding rocky objects originated from BP Psc's giant star disk material, then 
they likely formed through similar pathways as planetesimals forming around young stars. 
Another comm on indicator of pla netesimal formation in young stars is grain growth to mm- 
cm sizes (e.g., iNatta et all 120071 . and references therein). Is there any evidence for grain 
growth in the disk around BP Psc? Z08 report 880 /im measurements of BP Psc and 
interpret the continuum flux as originating from a population of cooler or potentially larger 
grains. The only other longer wavelength data point that exists for BP Psc that could be 
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used to constrain grain growth is an upper limit reported at 1 mm (Z08). This upper limit, 
however, does not allow one to place any stringent limits on the maximum grain size within 
BP Psc's dusty disk. Limits on the growth of grains in the disk around BP Psc will require 
more constraining observations at wavelengths longer than 880 /im. 

There are some pit-falls in this proposed dust enrichment model. If the colliding plan- 
etesimals are left over from BP Psc's initial planetary system, then there must be a large 
population of trojan-like bodies whose orbits evolve to bring them closer to the giant planet. 
In theory, interactions between such rocky bodies and BP Psc's gaseous disk could alter the 
rocky bodies' orbits and bring a sufficient number of them close enough to the giant planet 
to experience enhanced collisional activity (the likes of which must produce ~1 Lunar mass 
of small crystalline grains; see Table 1). Of course, if the planetesimals are instead being 
formed in-situ in the g iant star disk, then this problem disappears. Another problem is that 
Bouwman et all (120031 ) do not provide a convincing argument as to why one should expect 
isolated regions of nearly 100% crystalline dust in the disk atmosphere when the colliding 
bodies that produce the dust are more compositionally diverse (see the middle panel of 
Figure 5, Figure 10, and Table 3 in lBouwman et q/.ll2003l ). 



6.2.2. Giant Planet Disk Shocking 



Another possibility is that the giant planet generates shocks in the disk that sufficiently 
increase the disk temperature to anneal amorphous grains. Many theoretical works exist that 
predict shocks in the edges of a disk gap as a result of tidal in t eractions between the (gap- 
open ing) giant planet and the disk (e.g., iTakeuchi et a/.lll996l ; iLubow et a/.lll999t iRafikov 
20021 ). These shocks could act to heat small, amorphous grains to sufficient temperature s 
such that thermal annealing into crystalline species would result (IHarker Deschl |2002| ). 
The most interesting aspect of this model is the dust mineralogy it predicts. Annealing 
Mg-rich amorphous grains results in the generation of forsterite and Si02 as a by-product, 



while reaction of forsterite and SiC>2 results in the production of enstatite (jBouwman et al. 



200 ll . and references therein). Thus, giant-planet induced shocks in the disk are capable 
of reproducing the observed small grain mineralogy. F urthermore, shocks are the favored 
production method for cristobalite (jSargent et a/.ll2009bl ). the most likely polymorph of SiC>2 
observed in BP Psc's disk. 

After shocks produce the crystalline material in BP Psc's disk, the processed dust must 
be relocated from the disk mid-plane into the disk atmosphere where it is observed. Assuming 
shock-generated crystalline material enters the disk gap, it may then be orbitally perturbed 
(both radially and vertically) by the giant planet such that it may be directly irradiated 
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by stellar light and driven out of the disk mid-plane by radiation pressure. The flared disk 
geometry will allow these liberated grains to be re-captured by the disk, although at larger 
separations than where they originated. Thus, the disk atmosphere will become polluted 
with highly crystalline material as is observed for BP Psc. 

This is our preferred model as it naturally accounts for both BP Psc's disk structure 
and mineralogy. 



7. Conclusions 

We have obtained and analyzed Spitzer IRS spectra for the enigmatic dust and gas- 
enshrouded star, BP Psc. Disk structure analysis suggests that BP Psc's disk is unlike any 
classical T Tauri star disk, but is instead similar to disks around Herbig Ae/Be stars. Such 
a result suggests that BP Psc is itself as luminous as these pre-main sequence, intermediate 
mass stars. The suite of evidence presented in Z08 requires that BP Psc — if it is indeed 
as luminous as a typical HAeBe star — be a first-ascent giant star. The match of BP 
Psc's disk structure with that of the HAeBe star HD 100546 suggests that BP Psc hosts a 
gapped disk. Evidence for a giant planet as the generator of the gap in HD 100546 's disk 
suggests that BP Psc could also have a giant planet orbiting within its disk. Dust model 
fits to the IRS spectrum show that BP Psc's flared disk atmosphere contains chemically and 
thermally processed, cool, small dust grains. Our preferred model to explain the location 
and mineralogy of these grains relies on a giant planet orbiting BP Psc. Disk shocks that 
result from disk-planet interaction generate the highly crystalline dust which is subsequently 
removed from the disk mid-plane and relocated into the di sk atmosphere where it is observed. 



We briefly note, as suggested by iBouwman et all ( 120031 ). that this highly crystalline dust 



could be the seed material that forms cometary bodies like Hale-Bopp. If this is the case, 
and if BP Psc is indeed a first-ascent giant star, then we could be witnessing the initial 
formation stages of a second generation of cometary bodies around BP Psc. 
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Fig. 1. — All Panels: Reduced, extracted, and calibrated IRS spectrum of BP Psc. The 
spectrum S/N per pixel is >100. Top Panel: The ~5-12 jj,m wavelength range. Weak 
features in the ~7-8 jum range are dubious as this is an order-overlap region. Note, however, 
that the ~6.5 /im absorption/emission complex is real and is confirmed in the spectra of 
both nods. Middle Panel: The ~12-18 jj,m. wavelength range. Weak features in the ~14-15.2 
/im range are dubious as this is an order-overlap region. Bottom Panel: The ~ 18-35 jum 
wavelength range. Weak features in the ~ 19-20 range are dubious as this is an order- 
overlap region. Data beyond 35 //m are excluded due to the numerous bad pixels that reside 
in that region of the LL chip. 
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Fig. 2. — Top Left Panel: Spectral energy distribution of BP Psc (from IZuckerman et al. 
20081 ) with the reduced, extracted, and calibrated IRS spectrum and MIPS 24 /mi mea- 
surement overlaid. Top Right Panel: S pectral energy dist ribution of HD 100546 with the 
calibrat ed ISO-SWS spectrum overlaid (IMalfait et a/.lll998l ). A 10000 K stellar photosphere 
model (IHauschildt et a/.lll999[ ) has been fit to the BVRI data points. Bottom Left Panel: 
BP Psc measurements after subtracting the 5000 K photospheric model. The residual disk 
flux is fit with two blackbodies (brown, dashed lines) having temperatures of 1500 K and 
200 K, respectively. The 200 K blackbody displayed is similar to the ~200 K blackbodies 
used in fitting the dust continuum emission in the IRS spectrum (see Figure E]). Bottom 
Right Panel: HD 100546 measurements after subtracting the 10000 K photospheric model. 
The residual disk flux is fit with two blackbodies (brown, dashed lines) having temperatures 
of 1500 K and 200 K, respectively. 
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Fig. 3. — IRS spectrum of BP Psc; the spectrum S/N per pixel is >100. Model fits to 
the data (upper black curve) are overplotted with individual dust species' emissivities offset 
by negative flux for clarity. From bottom to top the models are for: cristobalite (black 
solid curve), pyroxene (purple dashed curve), enstatite (green dashed curve), forsterite (blue 
dashed curve), blackbodies (three brown dotted curves) having temperatures of 175, 250, and 
1000 K respectively, and the sum of all the aforementioned curves (red dot-dashed curve). 
See Table [1] for the dust parameters. 



1.0 



6 



8 10 

Wavelength (microns) 



12 




14 15 16 17 18 

Wavelength (microns) 




18 19 20 21 22 23 

Wavelength (microns) 



Fig. 4. — Regions of BP Psc's IRS spectrum having features from SiC^; the spectrum S/N 
per pixel is >100. Overplotted are final dust models using cristobalite (green, dashed line; 
see also Figure [3]), silica (red, dash-dot line — see Section 1X21) . and a model with no silica 
or cristobalite (gold, dash- dot- dot- dot line). Fits using cristobalite result in lower \ 2 values 
than those employing silica and no silica or cristobalite. Cristobalite performs better around 
12.7 jum, 16 jUin, and 20 fim. 
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Table 1. Dust Component Parameters 


Dust Component 


Stoichiometry 


Temperature 
(K) 


Mass Fraction a 
(%) 


Mass b 
(Lunar Masses) 


Forsterite 
Enstatite 
Olivine 
Pyroxene 
Cristobalite 
Continuum 1 
Continuum 2 
Continuum 3 


Mg 2 Si0 4 c 
MgSi0 3 c 
Mg 2 Si0 4 d 
MgSi0 3 d 
Si0 2 c 


75±25 
75±25 
150±25 
150±25 
250±50 
10001™ 
25011 
175+ 25 


59.41^ 

36.2+l| 
06+ 035 

u - uu -0.06 

4-o±8J 

0.03j£°» 


1.2 
0.6 

e 

6xl0~ 2 
6xl0" 4 

600 f 



a Fraction of mass in 0.1-1.0 //m, solid-state emitting grains detected in our IRS 
spectrum only. The reported percentage does not include the mass contained in the 
material responsible for the continuum emission. 

b Assuming the mass of the moon is 7.35 xlO 25 g and that BP Psc is a giant star 300 
pc distant from the Earth (Z08). For an assumed distance of 80 pc, appropriate if BP 
Psc were a T Tauri star, the mass per grain species is reduced by a factor of ~15. 

c Crystalline grains (see Section |3~2|) . 

d Amorphous grains (see Section [372]) . 

e We do not consider amorphous olivine as significantly detected in our final x 2 model 
fit. 

f This mass is from Z08 where they computed the disk dust mass assuming half of 
the 880 /xm detected flux towards BP Psc comes from a region of the disk having 
Tdust= 36 K (fits to our I RS spectrum do no t require that an additional ~36 K dust 



region exists; Figure 4 of iMalfait et a/.lll998l shows how the 33 /iin forsterite feature 



is sensitive to cool material). If there is no 36 K dust, then the mass listed for the 
continuum component could be a few times smaller. 



